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Abstract 

Based on the Boltzmann-Uehling-Uhlenbeck transport model coupled with the Lorentz 
force equation, we studied nucleus-nucleus collision in strong magnetic field. We find 
that neutrons and protons can be separated from a nucleus by strong magnetic field 
and neutron-rich high density nuclear matter and low density proton collectivity 
matter can be formed during nucleus-nucleus collision. The electric field produced 
by proton collectivity can accelerate proton and charged meson up to very high en- 
ergies. Besides the studies of isospin physics such as symmetry energy, these results 
may help us to understand the acceleration mechanisms of high energy charged 
particles in the cosmic rays. 
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The question of the origin of cosmic rays continues to be regarded as an 
unsolved problem even after almost one century years of research since the 
announcement of their discovery in 1912 [1]. The cosmic rays of extremely 
high-energy pose a serious challenge for conventional theories of origin of cos- 
mic rays based on acceleration of charged particles in powerful astrophysical 
objects. The question of origin of these extremely high-energy cosmic rays is 
currently a subject of much intense debate and discussions [2]. On the origin 
of energetic cosmic-rays in the universe, except through decay of sufficiently 
massive particles originating from processes in the early universe, there are 
basically two kinds of acceleration mechanisms considered in connection with 
cosmic rays acceleration, i.e., direct acceleration of charged particles by an 
electric field and statistical acceleration in a magnetized plasma. In the direct 
acceleration mechanism, the electric field in question is generally due to a ro- 
tating magnetic neutron star (pulsar) or, a (rotating) accretion disk threaded 
by magnetic fields, etc. pQ. In this article, we show that nuclear collision in 
strong magnetic field can produce proton collectivity in space, protons or other 
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charged particles thus can be accelerated directly by electric fields, which are 
produced by positive charges, not by changing magnetic fields. 

The condition of strong magnetic field may exist in the universe, such as 
white dwarfs, neutron stars, and accretion disks around black holes, and the 
maximum value of magnetic fields in the universe may reach 10 20 ~ 10 42 G [3J. 
And with the rapid development of laser technology, obtaining strong magnetic 
field artificially in terrestrial laboratory also may be possible [UEj. Although 
nuclear collisions exist extensively in the universe, up to now, nuclear reactions 
in the strong magnetic field were seldom reported. The separation of neutrons 
and protons from a nucleus while nucleus-nucleus collision in strong magnetic 
field may shed light on the acceleration mechanisms of high energy charged 
particles in the cosmic rays. 

Nuclear matter studies can be carried out semiclassically following the so 
called nuclear pasta phases [6fTf8f9] . In this study we also use the semiclassi- 
cal isospin-dependent Boltzmann-Uehling-Uhlenbeck (BUU) transport model 
[TDPTj . which is quite successful in describing dynamical evolution of nuclear 
collision. The BUU equation describes time evolution of the single particle 
phase space distribution function f(r,p,t), it reads: 
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f(r,p,t) can be viewed semi-classically as the probability of finding a particle 
at time t with momentum p at position r. The mean-field potential U depends 
on position and momentum of the particle and is computed self-consistently us- 
ing the distribution functions f(r,p, t). The collision item Illusion on the right- 
hand side of Eq. (PQ) governs the modifications of f(r, p, t) by elastic and inelas- 
tic two body collisions caused by short-range residual interactions [T2]. The 
proton and neutron densities of colliding nuclei are given by Skyrme- Hartree- 
Fock with Skyrme M* force parameters [13]. The isospin dependence is in- 
cluded in the dynamics through nucleon-nucleon collisions by using isospin- 
dependent cross sections and Pauli blocking factors [Hj. The isospin and 
momentum-dependent mean field potential used is [15] 
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In the above equation, 5 = (p n — p p ) / (p n + p p ) is the isospin asymmetry 
parameter, p = p n + p p is the baryon density and p n , p p are the neutron and 
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proton densities, respectively, r = l/2(— 1/2) for neutron (proton) and r ^ t', 
a = 4/3, / r (r, p) is the phase-space distribution function at coordinate r and 
momentum p. The parameters A u (x), Ai(x), B,C T)T , C T y and A were set by 
reproducing the momentum-dependent potential U(p, 5, p, r) predicted by the 
Gogny Hartree-Fock and/or the Brueckner-Hartree-Fock calculations [T6"fT7] . 
the saturation properties of symmetric nuclear matter and the symmetry en- 
ergy of about 31.6 MeV at normal nuclear matter density po = 0.16 fm~ 3 . 
The incompressibility of symmetric nuclear matter at normal density is set 
to be 211 MeV. The parameters A u (x) and A\(x) depend on the parameter x 
according to 

2B IB 

A u (x) = -95.98 x, Ai(x) = -120.57 + x, (3) 

a + 1 o + 1 

where B = 106.35 MeV. A = p° F is the nucleon Fermi momentum in symmetric 
nuclear matter, C TjT > = —103.4 MeV and C r r = —11.7 MeV. The C T y and C T)T 
items are the momentum-dependent interactions of a nucleon with unlike and 
like nucleons in the surrounding nuclear matter. The parameter x is introduced 
to mimic various density-dependent symmetry energies E sym (p) predicted by 
microscopic and phenomenological many-body approaches. Because we do not 
study nuclear symmetry energy here, in this article we just let the variable x 
be 1 [18]. The isoscalar part {U n + U p )/2 of the single nucleon potential was 
shown to be in good agreement with that of the variational many-body calcu- 
lations and the results of the BHF approach including three-body forces. The 
isovector part (U n — U P )/2S is consistent with the experimental Lane potential 
[H] . We use the isospin-dependent in-medium reduced nucleon- nucleon elastic 
scattering cross section from the scaling model according to nucleon effective 
mass. For in-medium nucleon-nucleon inelastic scattering cross section, we at 
present use the free nucleon-nucleon inelastic scattering cross section [T9j . 

Updates of nucleonic momentum are generally owing to momentum and spatial 
location dependence of its mean-field potential U, decided by the gradient 
force VU in Eq. ([1]). Besides the gradient and Coulomb forces added on the 
charged particles, momentum of charged particle also changes owing to the 
Lorentz force. For the additional magnetic filed force of charged hadron, we 
use the Lorentz force equation 

~f = q lf x if. (4) 

Where if is the velocity of charged particle and if is the additional magnetic 
field. In the practical calculations we add the Lorentz force to the Coulomb 
and mean-field gradient forces. Because the strength of magnetic field in this 
study is huge, we in the present simulations use the relativistic form of the 
Lorentz force equation and let the time step interval dt of updates of particle's 
phase space information be a very small value (dt = 0.0025/m/c). 
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Studying nucleus-nucleus collision in strong magnetic field, we should firstly 
get a picture of how the collision evolves with or without magnetic field. Fig. 1 
shows nuclear collision process with and without strong magnetic field. From 
upper panels we can see that, without magnetic field, protons and neutrons 
have almost the same mode of motions. Protons and neutrons go through the 
same compressions and inflations. These are normal knowledge of nuclear re- 
action. With strong magnetic field, however, the whole situation is changed. 
From the lower panels we can first see that, protons in the target and projectile 
refuse to collide and look like keeping still collectively owing to the Lorentz 
forces added on the protons. Whereas the neutrons in the target and projectile 
trend to collide, just like the case without magnetic field. The colliding neu- 
trons form high density neutron matter transiently. The "still" protons of the 
target and projectile form low density proton matter. We dub neutron matter 
or proton matter asymmetric nuclear matter. If the strength of magnetic field 
is smaller than ~ 10 17 Tesla, the Lorentz force can not overcome the mean- 
field gradient force. The separation thus can not happen. With such strong 
magnetic field, the protons almost keep still in space, i.e., fixed in coordinate 
space by the magnetic field and can be kept macroscopically. The separation 
of neutrons and protons from a nucleus with strong magnetic field may shed 
light on the possible origin of cosmic rays. If the magnetic fields are not ho- 
mogeneous, such as the wandering magnetic fields, the separation of neutrons 
and protons within a nucleus can be kept ultimately. And the separating neu- 
trons and protons move respectively. In the universe, the collision between 
nuclei may be replaced by physical collisions between stars [2"(J(21f22] . Colli- 
sions among light nuclei in stars or interstellar matter with magnetic field are 
similar to the 124 Sn+ 112 Sn collision here. In case the magnetic field decreases 
or disappears, the existing proton collectivity disperses promptly owing to the 
Coulomb actions. 

To show the acceleration mechanisms of high energy charged particles in the 
cosmic rays by the electric field, we draw the sketch-map of Fig. 2. The left- 
hand side of Fig. 2 is the sketch-map of proton collectivity forms in magnetic 
field via large number of "nucleus-nucleus collision" (in fact it is not nucleus- 
nucleus collision, but mainly is neutron-neutron collisions and protons keep 
"still"). When two dense stars with strong magnetic field collide, as shown 
in Fig. 2, large number of proton collectivity can be formed. They are fixed 
in certain space by the strong magnetic field so that the coulomb potentials 
are stored. While without strong magnetic field, however, this configuration 
changes promptly. The right-hand side shows the case without strong magnetic 
field. Protons in the collectivity disperse in all directions with the actions of 
the Coulomb potential. The first escaped protons or the boundary protons will 
be greatly accelerated, so they possess large velocities. The inner protons will 
have small velocities. Some escaped protons can be accelerated again and again 
here and there. This may be a reason why the cosmic-ray composition (above 
1.6 EeV) is proton-dominated [2U|24f25j . The energy of proton accelerated can 



4 



be roughly calculated via 

Here U is the Coulomb energy (in unit of GeV), Q (in unit of elementary 
charge) and R (in unit of fermi) are the accelerating charges and mean distance 
between accelerated proton and accelerating charges, respectively. If the value 
Q is large enough, the proton then can obtain a energy of several or even 
tens/hundreds of GeV. Assuming the accelerated proton mean R is 100 fm 
and proton density in collectivity (for proton collectivity, we suppose a semi- 
spherical shell distribution) is the same as that in 208 Pb, the proton then can 
obtain a maximal energy of about 0.03 x R/ fm x GeV = 3GeV. Practically, 
scale of proton collectivity formed in strong magnetic field may be much larger 
than 100 fm, the proton energy is thus also much larger than 3GeV. Note here 
that proton may be accelerated many times as shown on the right-hand side of 
Fig. ?? owing to the winding magnetic field. Therefore the energetic cosmic-ray 
protons or other charged particles may be accelerated by electric field. As for 
other cosmic-ray nuclei, they can be accelerated directly by the electric field, 
having no colliding separation of neutrons and protons via nucleus-nucleus 
collision. 



Because of the separation of neutrons and protons within a nucleus while 
nucleus-nucleus collision in strong magnetic field, the compressed nuclear mat- 
ter must be neutron-rich matter. Fig. 3 shows neutron to proton ratio n/p of 
such high density nuclear matter formed during the collision with strong mag- 
netic field. Here the n/p of high density nuclear matter formed can reach 4 (if 
the colliding nuclei are far apart, the n/p of compressed nuclear matter can be 
even more large), is much larger than obtained from general nuclear reaction 
in terrestrial laboratory, in which the n/p of transiently formed nuclear matter 
is about 1. If in the universe, somewhere there are strong magnetic fields as 
large as 10 17 Tesla, one then can possibly in heaven use such asymmetric nu- 
clear matter to study isospin physics (physics relevant to unequal numbers of 
neutrons and protons), such as nuclear matter symmetry energy (energy rel- 
evant to the changing of mean energy per nucleon owing to unequal numbers 
of protons and neutrons, is a main subject of isospin physics), which is crucial 
for understanding many interesting issues in both nuclear physics and astro- 
physics p6l27]l2 8 29 30 3~Tf32] and has been regarded as the most uncertain 
property of dense neutron-rich nuclear matter [33|3l] . The advantage of using 
such asymmetric nuclear matter to study isospin physics is that such asym- 
metric nuclear matter formed in strong magnetic field via nuclear collision has 
huge asymmetry, effects of isospin are thus enlarged remarkably. Ever since a 
long time ago, the small asymmetry of compressed nuclear matter obtained 
in terrestrial laboratory has been one of the most troublesome factors (which 
is owing to small resolution of isospin effect caused by small asymmetry of 
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nuclear matter produced through ordinary nuclear collision) to study isospin 
physics, especially for nuclear symmetry energy. 



Pion production in nuclear reaction nowadays has attracted much attention 
in nuclear physics community [28,29,19]. One important reason is that pion 
production is connected with the high density behavior of nuclear symmetry 
energy [53] • Here we only study pion meson in original cosmic rays. Fig. 4 shows 
evolution of n~ /ix + ratio with and without magnetic field, using 124 Sn+ 112 Sn 
collision at 400 MeV/nucleon as a example. At early stage of nucleus- nucleus 
collision in the magnetic field, as shown in Fig. ??, we can see that neutron- 
neutron collisions can happen while protons nearly do not take part in the 
collisions. We thus see a very large tt~/tt + ratio with magnetic field (7r~'s are 
mainly from neutron- neutron collisions, 7r +, s are mainly from proton-proton 
collisions, neutron-proton collisions produce roughly equal numbers of ir~ and 
7r + ). In the universe, with strong magnetic field, such behavior can happen dur- 
ing nucleus-nucleus collision. Maybe some of n mesons, especially ir~ mesons 
in the cosmic rays come from such nucleus-nucleus collision, and they 
celerated (the acceleration mechanisms may be that pions are attracted by 
proton collectivity and then their trajectories are bent by magnetic field) by 
the electric field produced by proton collectivity originating from nucleus- 
nucleus collision in strong magnetic field. 

In conclusion, we do a study of nuclear collision with strong magnetic field. 
It is found that nuclear collision in strong magnetic field (roughly ~ 10 17 
Tesla) causes the separation of neutrons and protons from a nucleus. The 
formed proton collectivity can produce electric field and accelerates proton 
or charged pion meson up to very high speeds. Nuclear collision in strong 
magnetic field can also produce asymmetric nuclear matter transiently. The 
formed asymmetric nuclear matter may be used to study isospin physics such 
as nuclear symmetry energy someday in both heaven and terrestrial laboratory. 
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